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for several values of the latter, and the vacuum energy and the chiral condensate 



[ potentials is addressed and it is demonstrated that the number of replicas is infinite 

Q^' for any a, as a consequence of the peculiar behaviour of the quark self-energy in the 

. infrared domain. 

: 
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I. INTRODUCTION 

The problem of the spontaneous breaking of chiral symmetry (SBCS) and its relation to 
confinement is one of the QCD cornerstones. Although the basic ideas of SBCS are already 
a textbook topic, this problem still lies at the crossroad of many studies and approaches to 
QCD. In this paper we exploit the potential model for QCD, whose origins can be traced back 
to QCD in the truncated Coulomb gauge and which is proven to be successful in studies 
of the low-energy phenomena in QCD (see, for example, This class of models can 

be indicated as Nambu-Jona-Lasinio-(NJL-)type models [2] with the current-current quark 
interaction and the corresponding formfactor coming from the bilocal gluonic correlator. 
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The standard approximation in such type of models is to neglect the retardation and to 
approximate the gluonic correlator by a confining potential of a certain form. The power- 
like potentials, which are the most natural candidates for the role of the confining force, are 
the subject of the present investigation. In the course of this paper we 

• revisit the problem of SBCS for the for power-like confining potentials V{r) = i^'g^^r" 
with a ^ 0; 

• restrict the range of allowed a's and, for several values of the latter, find numeri- 
cally solutions to the corresponding mass-gap equation, as well as the vacuum energy 
density, and the chiral condensate for the chirally noninvariant vacuum of the theory; 

• study in detail the problem of existence of the replica solutions to the mass-gap equa- 
tion for power-like confining potentials. 

The problem of instability of the chirally invariant vacuum for power-like confining poten- 
tials was studied in detail in the middle of 80's by the Orsay group and such an instability 
was proved for the range ^ a < 3. For numerical studies the harmonic oscillator type 
potential, a = 2, was chosen by these authors, as well as by the Lisbon group |^], and a 
set of results for hadronic properties were obtained in the framework of the given model. 
In this paper we study the mass-gap equation for an arbitrary value of a ranging from 
to 2, with special attention payed to a's close to unity, since the linearly rising potential is 
known to be preferred by phenomenology as the most successful candidate for the confining 
force giving the correct Regge trajectories behaviour, possessing a clear connection to the 
QCD string, and so on (see p] and references therein). It is also claimed to be singled out 
by the lattice calculations. We exclude the region a > 2 since the corresponding mass-gap 
equation diverges for such a's. On the other hand, it would be hard to justify the use of 
such a strong confining force in phenomenological models for QCD. We solve the mass-gap 
equation explicitly for several values of a from the allowed region and demonstrate that the 
chiral angle, the vacuum energy density, and the chiral condensate are smooth slow func- 
tions of the form of the confining potential, so that the results obtained for the potential 
of a given form — the linear confinement being the most justified and phenomenologically 
successful choice — have a universal nature for any quark-quark kernel of such a type. 

Following the set of recent publications devoted to possible multiple solutions for the 
chirally noninvariant vacuum in QCD j^, Q| (see also where a similar conclusion was made 
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in a different approach), we address the question of rephcas existence for various power laws 
r". We find that for a whole range of allowed powers, ^ a ^ 2, replica solutions do exist, 
similarly to the case of a = 2 studied in detail in P, 0]. We give the profiles of several 
replicas for the linear confinement and argue that the number of such solutions is infinite 
for any power a, including the weakest, logarithmic, potential which corresponds to a = 0. 
We argue that the source of replicas is the infrared behaviour of the single-quark self-energy 
— the dressed quark dispersive law E{p) — which, for small values of the quark momentum 
p, becomes a negative sharp function of p thus enabling fast oscillations of the chiral angle 
with the frequency increasing with vanishing momentum p. Since this property of the quark 
dispersive law is expected to be an integral part of any confining interaction, then we confirm 
the conclusion made in that ". . . across all these different quark kernels, the existence 
of vacuum replicas should constitute the rule rather than the exception." We argue that, in 
real QCD, with the confining interaction flattening at large distances due to the effect of 
the QCD string breaking, the number of replicas becomes finite. We find the parameter of 
the SBCS given by the replica solutions to decrease fast with the number of nodes of the 
chiral angle so that one has a well defined perturbative series in replicas and, therefore, the 
account only for the first replica may be sufficient in many phenomenological applications 
(the details of the formalism which allows one to incorporate replicas into quark models can 
be found in ref. |3]). 

The paper is organised as follows: In the second section we give the necessary details of 
the formalism and derive the mass-gap equation for the power-like confining potential, which 
is studied in detail in the third section, first, qualitatively, then quantitatively, and, finally, 
numerically. The mass-gap equation is solved numerically for several values of a and the 
chiral condensate and the excess of the vacuum energy density over the trivial solution are 
calculated for the found solutions. In the fourth section, devoted to replicas, we demonstrate 
how an infinite number of solutions to the mass-gap equation appears and build explicitly 
a couple of replicas for the linear confinement. Our conclusions are the subject of the last, 
fifth, section. 
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II. THE MASS-GAP EQUATION 



The chiral model which we use for our studies is given by the Hamiltonian with the 
current-current interaction parametrised by the bilocal correlator K^j^^, 



H 



d'x^{x,t) (-^7- v)^(f,t) + ^ / d'xd'yj;{x,t)K;l{x-y)J'Mt), (1) 



where the quark current is J^{x,t) = 'ip(x,t)'j^^ip{x,t), and the gluonic correlator is ap- 
proximated by a potential, 



Ktix - y) = g,o9uoS^%{\x - y\) 



(2) 



with 



Voi\x\)=Kr'\xr. (3) 

In order to include the logarithmic potential into consideration the obvious modification 
of the potential is needed: 

{Ko\x\r-i 



a 



Ko\n{Ko\x\). 



(4) 



The model contains the only dimensional parameter — the strength of the confining_force 
Kq. For further convenience we shall consider a modified version of the potential (jSj) |2|, 



(5) 



where m plays the role of the regulator for the infrared behaviour of the interaction. The 
limit m — > is understood. 

The standard technique used in such models is the Bogoliubov-Valatin transformation 
from bare to dressed quarks parametrised by the chiral angle — the main entity defining 
the chiral symmetry breaking, the structure of the BCS vacuum of the theory, as well as 
the properties of the hadronic states built over this vacuum For the application of 

this technique to the two-dimensional QCD see the papers 10 1. We choose the following 
parameterisation: 

I I 

^(P) = 775 



(6) 



V2 
1 
V2 



1 + sin ifip + 1 - sin (pp (ap) 
^\ + sini/jp - - siny^p (ap) 



m(0). 



(7) 
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b^ip, t) = e'^^^'b^ip, 0), t) = e'^'^'d^i-p, 0), (8) 

where Ep (the shorthand notation for E{p)) stands for the dispersive law of the dressed 
quarks, the chiral angle (p{p) (we also use the shorthand notation ipp for it) varies in the 
range — 7r/2 < (/9p ^ 7r/2 with the boundary conditions v^(0) = 7r/2, ip{p oo) 0. 

The Hamiltonian normally arranged in the basis (jH)) splits into the vacuum energy, 
the quadratic and the quartic parts in terms of the quark creation/annihilation operators. 
For the vacuum energy density one has 

^vacM = ^{Q\TH[ip]\0) = -f y -0^{A{p)smipp+[B{p) +p]cos<^p), (9) 

where V is the three-dimensional volume; the degeneracy factor g counts the number of 
independent quark degrees of freedom, 

g = {2s + l)NcNf, (10) 

with s = i being the quark spin; the number of colours, Nc, is put to three, and the number 
of light flavours, Nf, is two. Thus we find that g = 12. The auxiliary functions A{p) and 
B{p) are defined as 

A{p) = -Cf j j^Mp- ^) si^'/'^' (11) 

B{p) =p + -Cp j ^ (p ■ k)V,{p- k) cosy.,, (12) 

where Cp = | is the 5'f/(3)c Casimir operator in the fundamental representation. The 
actual form of the chiral angle is such that the quadratic part of the normally ordered 
Hamiltonian diagonalises, or alternatively, the vacuum energy takes its minimal value. The 
corresponding equation 

= 0, (13) 

known as the mass-gap equation, reads 

A{j)) cos ifp — B{p) sin ipp = 0. (14) 



For the generalised power-like potential one can find 

1"^ dkkcos[{a + l 



K^^^ /•°° dA;A;cos[(a + l)arctan^] ^ , 
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B{p) =p-CFT{a 



K^+^ r dk 



271 



apk cos[(a + 1) arctan cos[(q; — 1) arctan 



[m^ + (A; — pY] 2 



(k — p) sin[a; arctan 



cosv^fc, 



{a — l)[m? + {k — pY] 2 

(16) 



[m? + {k — 2^ 

where, for the sake of convenience, we continued the integral to the negative values of 
k assuming cos = —cos sin^j.^ = siny^^ (the most natural realisation of these 
conditions can be achieved in terms of some even function m„, such that sino9„ = —7^S== 



cosyjp 



which plays the role of the effective mass of the quark). Consequently, 



the mass-gap equation takes the form: 



p'sm^, = CFKt^T{a) / — 



dk J apk cos [(a + 1) arctan ^— 



271 



[m? + [k — pY] "2 



+ 



cos[(a — 1) arctan [k — p) sin[a arctan 



{a ~ l)[m^ + {k - pf]"^ [w? + {k-py]2 

and it is the main object of our studies. 



sin[99p - ifk] 



cos ifk sm yjp 



(17) 



III. INVESTIGATION OF THE GENERAL FORMULA 



A. Qualitative analysis 

As the first step in studies of the general formula (fTTj) . we perform its simple qualitative 
analysis. Using the techniques described in j^, we assume that a solution ^^^\p) to this 
equation exists, and the vacuum energy is minimal on this solution. If the function ip'^^^p/A) 
with an arbitrary scale parameter ^ A < 00 is substituted into the vacuum energy Q, 
then the function £^s.c{^) must reveal a minimum for A = 1. Moreover, the corresponding 
minimum should lie lower than the one for the trivial solution ^{p) = 0, when there is 
no dressing of quarks and the chiral symmetry is unbroken. If the regulator m is removed 
from the vacuum energy functional, then the strength of the potential remains the only 
dimensional parameter in the theory, so that, after a proper rescaling of the integration 
variables, one arrives at the simple formula, for an arbitrary number of spatial dimensions 
d, 

= C.A^^' + C^K^+^A"-^, (18) 
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where Ci and C2 are two constants independent of A which are interrelated by the constraint 
dEyii^{A) / d A\A=i = 0. The first term in Eq. p8|) comes from the kinetic energy, the second 
term is due to the interaction. The following four situations are possible: i) < a < d, ii) 
a > d, and two boundary cases, iii) a = d, and iv) a = 0. In the first case the vacuum 
energy has a double- well form with two minima: trivial for A = and nontrivial for A = 1. 
The difference S^^dA = 1) — £^vac(^ = 0) is negative, so that the chirally nonsymmetric 
nontrivial solution is energetically preferable. For the second case one has the interaction 
term in ()18p containing negative powers of A and, as a result, the trivial solution, with 
unbroken chiral symmetry and which corresponds to A = 0, possesses an infinite energy and 
therefore it does not exist. In the meantime, a nontrivial solution with A = 1 still may be 
present. The boundary case of a = d leads to the logarithmic dependence of the vacuum 
energy on the parameter A, 

^vac(A) = C,A''+' + C^K^^' In A, (19) 

so that qualitatively the same conclusion holds — the theory possesses only chirally non- 
symmetric phase. The two-dimensional QCD [9j is an example of the theory with such a 
logarithmic dependence (see ^] for the details). 

Finally, for the case iv), that is, for the logarithmic potential (@]) one has 

£..M) = CiA^+i + C2K0A'' In A, (20) 

where the logarithmic growth of the energy, when approaching the trivial solution limit A = 
0, is cancelled by the power factor A'^, so that both, chirally symmetric and nonsymmetric, 
solutions coexist in this case, similarly to other potentials with < a < d. 

Thus we conclude that once the power of the potential reaches the critical value equal 
to the number of spatial dimensions the behaviour of the theory changes drastically, the 
chirally symmetric phase being swept off. Meanwhile the qualitative analysis performed 
above ignored the problem of convergence of the integrals in the expression for the vacuum 
energy and in the corresponding mass-gap equation. It also fails to answer the question 
as to how many solutions to the mass-gap equation exist. In what follows we turn to the 
quantitative and numerical analysis of these problems. 
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B. Quantitative analysis 

Now we turn to the detailed analysis of the mass-gap equation (fTTj) . such as the problem 
of the convergence, the allowed region for a's, the dependence on the regulator m, and so 
on. 

First of all, one can easily check that the case a = brings no difficulties — the r.h.s. 
vanishes if the limit a ^ is taken naively, whereas to arrive at the mass-gap equation for 
the logarithmic potential one is to divide the r.h.s. by a (see Eq. (@])) which leads to a finite 
result after taking the limit a — ^ 0. 

For a = 1 the divergent term proportional to l/(a — 1) vanishes on the r.h.s. of Eq. p7|l 
since the cosine of the chiral angle is odd. An accurate expansion of this term for a ^ 1 
brings about logarithmic terms. 

Now let us check the largest value of a which does not lead to divergences in the mass- 
gap equation. When the regulator m tends to zero, the first term in the curly brackets in 
Eq. ()17|1 . formally, is the most singular term for ~ p, and it can be written as 

P rfA;pfccos[(a + l)arctan^] . . , 
In the region k ^ p the integrand admits an expansion in powers {k — p)": 



cos[(a + 1) arctan 
27r[m2 + {k- p^] 



™^ ip'(k-p)^'+ik-pY 



(22) 



where the first term in the curly brackets is odd and, therefore, it vanishes in the integral in 
k around k = p. The remaining expression, as well as other terms on the r.h.s. of Eq. (jl7|) . 
behaves as \k—p\^~'^ and the corresponding contribution to the mass-gap equation converges 
for a < 2, which is the upper limit of the range of valid powers a. For such a's, 

^ p -jj- 

arctan — > — sign(A; — p) + 0(m), (23) 

m m-+o 2 

and the regulator m can be removed from the mass-gap equation. For the case of a = 2 the 
entire r.h.s. of Eq. (fTTjl vanishes after the substitution ()23|) . To be more precise, one has to 
keep next terms in the expansion (j2Sl), with positive powers of the regulator m, reproducing 
the well-known representation of the delta function, 

6{p — k) = lim 



m- 



Ti w? + {k — pY ' 
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and its derivatives. We shall consider this special case separately. 

Notice that the proof of the chirally symmetric vacuum instability given in was based 
on the consideration of the self-energy functional = / -^^Voip — k){pk) oc which 
gives infrared divergent contribution to integrals in d^p if a; ^ 3. In the meantime, as demon- 
strated above, the requirement of finiteness of the mass-gap equation imposes a stronger 
restriction on a: a ^ 2. 

One encounters no more difficulties for a's within the interval ^ a < 2 and, when the 
regulator m is removed, arrives at the mass-gap equation in the ultimate form: 

dk (pksm[(fk 



p-^ sin = C^ir^+T(a + 1) sin 



27r 



\p — k 



a+l 



+ 



COS Sm yjp 

a — 1 



la-l 



A 



(24) 



\p — k\ 

We introduced the term with an arbitrary mass parameter A, in order to emphasise 
the convergence of the integral for a = 1. This extra term does not contribute to the integral 
due to the parity of cos (pk- 

In particular, the mass-gap equation for the logarithmic potential follows from plj) in 
the limit a — > 0, if the proper modification of the potential, given in Eq. (jH), is applied. 



p smyjp 



dk 



4 J-oo\p-k\ 



[pk sin[v9fc - V^p] - {p-kf cos v?fc sin yjp] 



(25) 



For the case of a = 1, that is, for the linear confinement, the formula well known in the 
U_..eadilyre,„duoed(.ee,fo.exa.p.,0), 



p svmpp 



CfKI 



CfKI 



dk 



dk 



pk 



, 27r \ {p- kf 
Ap^k"^ 



sin[v9fc - ipp] - 



sva.[(pk - Vp] - In 
2pk 



\p — k\ 



A 



COS(y9fc sm 



(26) 



{p + kf 



+ In 



p — k 



p + k 



cos (fk sm (pp 



2tT { (p2 - p)2 

To complete our investigation let us consider the case of the harmonic oscillator potential, 
a = 2. The Fourier transform of the potential V{r) = K^r"^ is the Laplacian of the three- 
dimensional delta function, so that the resulting mass-gap equation becomes differential 



a 



p^ sin Lpp = -CpKl [pVp + '^PV'p + sin 2^p\ . 



(27) 



Formally, Eqs. p7|) . ()24|) remain valid for —1 < a < 0, and the ultraviolet divergence is 
encountered for a = —1, that is, for the Coulomb potential. We disregard this region since 



10 



the resulting force fails to be confining. Thus the valid confining potentials, in momentum 
space, range as, symbohcally, {2'KfKQ5^^\p) < V{p) ^ {2TifKl/^5^^\p), that it, from 
constant to the harmonic oscillator potential, respectively. 

For numerical investigation of the mass-gap equation it is convenient to evaluate 
analytically the contribution of the strip \k — p\ < \ into the integral on the r.h.s., 

h = ^^^^^r(a + 1) sin ^[pV;' + 2p^; + sin2^,], (28) 

where this integral admits an extra contribution if the width of the strip is chosen different 
from A. For a = 2 the dependence of the strip integral ()28|) on A disappears and the full 
mass-gap equation (jTTj) is readily reproduced as a consequence of the delta functional form 
of the Fourier transform of the harmonic oscillator potential. 



C. Numerical analysis 

In this subsection we present the results of the numerical studies of the mass-gap equation 
(I24|l with ^ a ^ 2. The equation for the harmonic oscillator potential ()27|1 is studied 
in detail in the literature, so the interested reader can find the details, for example, in 
refs. ^, ^. In Fig. 1(a) we plot the profile of the nontrivial solutions ipo{p) to the mass- 
gap equation for several values of a. The solution for a = 2, found in refs. P, 0], 
is also depicted for the sake of completeness. In Fig. 1(b) we present the results for the 
chiral condensate (qq) = —-^ Io° P^sin(/?o(p) = — and for the excess of the vacuum 
energy A£^vac = ^^vac[<y5o] — i^vac[v^ = 0] = — over the trivial vacuum as functions of a. 
From Fig. 1(b) one can see that A^ ^ A^- in the whole range of allowed a's and that their 
dependence on the form of the potential is smooth, solutions for a around unity being 
quite close to one another (see Fig. 1(a)). Thus we conclude that the concrete form of the 
confining potential does not play a crucial role for the physics of chiral symmetry breaking, 
resulting only in minor numerical changes. In particular one can see that the qualitative 
behaviour of the solution is very stable against deviations of the potential from the purely 
linear form, usually adopted in phenomenological models and declared to be confirmed by 
lattice calculations. Therefore, at least SBCS and the spectrum of low-lying hadrons will 
not be strongly affected in case the behaviour of the confining potential is slightly changed, 



deviating from linearity, as suggested in [ll|. 
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(a) (b) 

FIG. 1: The solutions to the mass-gap equation (|24jl for a = 0.1, 0.3, 0.5, 0.7, 0.9, 1.0, 1.1, 1.3, 1.5, 
1.7, 1.9, and 2.0 (plot (a)); the curves localised closer to the origin correspond to smaller a's, and 
the mass parameters defining the chiral condensate, as (qq) = — A^, and the excess of the vacuum 
energy density over the trivial vacuum, for two quark flavours and three colours, as Afvac '- 
(plot (b)). All dimensional quantities are given in the units of Kq. 



IV. THE REPLICAS 



As argued in a sequence of recent papers jaLD], it is possible that ". . . the same ultraviolet 
behaviour, for instance for the quark propagator, bifurcates to different solutions when we go 
to the low-energy domain. . ." in QCD, and such a replica was discovered for a phenomenology 
inspired potential. Besides that, the whole infinite tower of excited solutions for the mass-gap 
equation (jTTj) for the harmonic oscillator potential was found in ref. and also confirmed in 
6|. With the general form of the mass-gap equation ()17p and ()24|) we are in the position to 
investigate the problem of replicas existence for various power-like confining potentials. We 
find that any power-like potential r", with the power ^ a ^ 2, maintains replicas. In the 
paper jsl a detailed analysis was performed for the harmonic oscillator potential mass-gap 
equation (^Tjl and the existence of an infinite tower of solutions was proved analytically. We 
failed to repeat this analysis for the general form of the mass-gap equation ()24j) since, in 
contrast to Eq. ()27|) . Eq. is integral with the coefficients tuned to provide the overall 
convergence of the integral, but leaving no hope to use any expansions under the integral. 
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Instead let us use an approximate method in order to demonstrate how rephca solutions 
occur for the mass-gap equation (j211)- To this end we use the parametrisation of the chiral 
angle through the effective quark mass and introduce a new function ipp-. 



sin Lpp = TUpDp, cos (fp = pDp, ipp = p sin ipp, = ^ p'^ + m^. (29) 

It is also convenient to use the following integral: 

]_^2.\k-p\^^^- 2r(« + l)sinf y_^^"^" ^^'^ • 

Then the mass-gap equation can be rewritten in a simple physically transparent 
form: 

[2E{p) + Cp/^o"+Vr] = 0, (31) 

where the operator E{p) has the meaning of the quark self-energy and it is given by the 
following expression: 



, 1 CfK^\, ^ . Tia r dk \p^Dp-kWk , pkDpD, 
E{p) = j, — — r(a + l)sm— / — 



• (32) 



Dp pWp ' ' 2 7_^27r 
The even function m(p) takes its maximal value at p = and then decreases rapidly. 
Thus, for momenta larger than some po ~ ??i(0), Eq. ()31|1 can be linearised by putting 
Dp^ ~ \p\ in the expression ((221) • Then the integral on the r.h.s. of Eq. fl32|) is easily 
evaluated and the self-energy takes the form: 

2C^-^o^^r(a + l)sinf 

Notice that for any confining potential, a ^ 0, the second term on the r.h.s. of the ex- 
pression dominates in the low-momentum region bringing a large negative contribution 
to the quark self-energy (the corresponding term becomes logarithmic, ~ KqIu^, for the 
potential (j3))). This general feature of the quark self-energy in the confining potential has 
been discussed in literature not once (see, for example, 0,^, or the papers |3| where the 
case of QCD2 is discussed in detail) and it is known to play a crucial role for the properties 
of the theory. Besides that, replicas exist in the theory also due to this feature of the quark 
self-energy E{p). In order to demonstrate this let us write the linearised mass- gap equation 
fl31|) in the form of a Schrodinger-type equation in momentum space. 



^1, , , , 4:C F K+^T( a) sin ^ 



CfK^'\x\'' + 2\p\ 



tt{2 - a)\p\' 



ipp = eipp, (34) 
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and notice that the hnearisation suppresses the first, positive, term on the r.h.s. of Eq. 
and enhances the second, negative, term. As a result, we are interested in the odd (see the 
definition (f^n|) ) eigenstates of the linear equation (jH^ with negative eigenvalues. Each such 
state indicates a solution of the full nonlinear mass-gap equation ()24j) and can be used as 
the starting anzatz for the iterative numerical method to solve the latter. Therefore, in 
order to find the number of solutions to the mass-gap equation (|^. it is sufficient to count 
the odd eigenstates, with negative eigenvalues, of the linear equation (jMj) . If the Bohr- 
Sommerfeld quantisation procedure is applied directly to Eq. (jHH), then the quasiclassical 
integral, Iwkb = f^"''''^ x(p)dp, diverges logarithmically at p = 0, where the chiral angle is 
not small anymore and the approximation D^^ ^ \p\ obviously fails. In this region, the 
sharp behaviour of the self-energy is smeared by the effective quark mass, /i ~ mp{p 0), 
which plays the role of the effective regulator, so that we modify Eq. (jH^ accordingly: 



ACpK^^ria) sin 



7ra 
2 



(35) 



7r(2-a)(p2 + ;x2)"/2 

The spectrum of eigenstates of Eq. (jH^j) starts, for = 0, at the bottom of the deep 
well described by the effective potential V{jp) = 2E{p), Eq ~ V{p = 0) 



and, for some rimax, reaches zero from below. Then the spectrum continues for positive 
eigenvalues up to infinity. Since the last negative eigenenergy ^nmax is small, then we expand 
the quasiclassical integral accordingly and find: 



Ko 



c^r(|)r(^) 



1 



a{l+a) 



ra)r(Tfe) 



vrn 



2"r {^) 



(36) 



where we expressed sin ^ through the Euler Gamma functions as 



• ^« a ir(|)r(i±^) 

sin — = 2"'^ , V / 

2 r(a)r(^) 



(37) 



Notice that in actuality the regulator /i itself is not a constant, but it rapidly decreases 
with n, and it is such that all eigenvalues of Eq. ()35|1 vanish for all n's. Then the corre- 
sponding eigenfunctions ipn provide, according to the definition (jSHI), the solutions to the 
exact mass-gap equation (j211). From Eq. one can easily see that 

l/a 



finia) = Kq exp (-C^vrn), Ca 



v^r(^) 
2"r {^) 



(38) 
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Using a more accurate expansion in Eq. ()36|) , one can find the correction 6^ to the leading 
regime (jHHj) . /in(tt) = -f^o exp (— CQ,(7rn + as it was done in P]. Unhke the leading 
logarithmic term, this correction is sensitive to the concrete form of the regularization of 
the potential V{p), and we do not give it here. 

From the formula ()H8|1 we conclude that, for any a, the mass-gap equation supports 
an infinite number of solutions which, in the low- momentum region {p <^ -^0)5 behave as 

^n{p ~ 0) = ^ - Ae^"(-+^") + . . . , (39) 

then oscillate n times and, finally, approach zero for the infinite momentum. The constant 
Ca is monotonously decreasing function of a in the whole interval ^ a ^ 2 and it varies 
from Co = I exp[|(7 — 1— ■?/'(|))] ^ 1.08 to C2 = ^ ~ 0.71. Numerically, the same behaviour 
of the solutions to the mass-gap equation (j24j) is observed — for smaller a's and higher 77,'s 
the chiral angle becomes steeper at the origin, in accordance with the obvious identification 
^'n\p=o = -1/yUn oc cxp (Cq[7™ + 5J) (see Figs. 1(a), 2). 

The formula (jnH|) is approximate since the dependence on fi in Eq. (jH^j) reproduces only 
the gross features of the quark self-energy, whereas the exact dependence of the self-energy 
on the chiral angle is, in turn, a consequence of the mass-gap equation, and the problem 
becomes self-consistent. To estimate the accuracy of the formula (|38|) we continue it to 
a = 2 and compare our result with the result found in ^|: 

C2 = ^ ^ 0.71, C2(ref. [3]) = ^ ^ 0.76, 
v2 v7 

that is, the error is about 6 ^ 7%. 

As an indirect confirmation of the conclusion made above let us mention that no critical 
phenomena are observed in the mass-gap equation for any a as well as in the limit a — 2, 
so that all properties of the mass-gap equation, including the infinite number of solutions, 
may be continued from a = 2 to smaller a's. 

In Fig. 2 we give the profiles of the ground-state, as well as of the first two replica solutions 
for the linear confinement, a = 1. Solutions with larger number of nodes are also available 
for numerical study, but the method has to be very precise since each new solution possesses 
oscillations squeezed to zero and seen only with the help of methods with a sufficiently high 
resolution. According to the formula pUj). the behaviour of these solutions at the origin can 
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FIG. 2: The three first solutions to the mass-gap equation H26j) corresponding to the ground-state 
BCS vacuum (n = 0) and to the first two replica states (n = 1,2). The momentum p is given in 
the units of Kq. 

be approximated as 



so that yU„_|_i(l)//i„(l) ^ exp (^) ~ 10, that is, the parameter of SBCS decreases about ten 
times for each next replica. 

In ref. Q a method is proposed which allows one to take the vacuum replicas into account 
in quark models. The generalisation of this approach to the case of many replicas is straight- 
forward and assumes the summation over contributions of all replicas. On the other hand, 
since the mass parameter of the solution, for example, /i„, decreases fast (approximately ten 
times for the linear confinement) for each next solution, then one has a sort of perturbative 
series with the well defined convergence parameter. Therefore it may be sufficient to consider 
only the first replica neglecting the contribution of higher replicas, starting from the second 
one, which are hard to distinguish numerically from the trivial solution ifp = 0. To have 
a good phenomenological description of quarkonia one should supply the purely confining 
potential by the short-range Coulomb interaction and, possibly, by a constant term in order 








(40) 
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to fit for the right value of the chiral condensate (an attampt to evaluate this constant from 
first principles was undertaken in Q]). Such a potential was considered in 0, and the 
solutions for the ground-state BCS vacuum as well as for one vacuum replica were found 
numerically. 



V. CONCLUSIONS 

In this paper we complete the study of the power-like potentials K^'^^r"' from the point 
of view of SBCS and the number of nontrivial solutions of the mass-gap equation. We use 
the constructive method to prove the chirally symmetric vacuum instability for such poten- 
tials, solving the mass-gap equation explicitly and calculating the vacuum energy for the 
corresponding solution. We establish the region of allowed powers of a for such confining 
potentials which lead to a convergent mass-gap equation and a finite excess of the vacuum 
energy density over the trivial solution with unbroken chiral symmetry. Thus, using simple 
qualitative arguments, we demonstrate that for ^ a < d, d being the number of spa- 
tial dimensions, at least two, chirally symmetric and nonsymmetric, solutions should exist, 
whereas for a ^ d the trivial solution possesses an infinite energy density and disappears. In 
the meantime, the restriction for the parameter a exists, ^ a ^ 2, which comes from the 
fact that the corresponding mass- gap equation should be convergent. We find numerically 
solutions to the mass-gap equation for various values of a from the allowed region and, for 
the found solutions, evaluate the vacuum energy density and the chiral condensate which 
are given by the same scale, ~ A^, and turn out to be slow functions of the parameter a. 

We address the question of existence of the second, third, and higher chirally nonsym- 
metric solutions of the mass-gap equation for power-like confining potentials and find that 
any potential with ^ a ^ 2 supports such solutions — the replicas — which is rather the 
rule for confining potentials and comes from the peculiar behaviour of the quark self-energy 
in the infrared domain. We find that the number of such replicas is infinite for any a and 
estimate the slope of the solutions to the mass-gap equation at the origin. 

Thus we dare predict the existence of replicas regardless of the explicit form of the 
confinement and of the details of the model used in calculations. In real QCD, with light 
quark fiavours and the quark-quark potential fiattening at large distances due to the effect of 
the QCD string breaking (such an effect can be taken into account in the interquark potential 
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through a coordinate dependence of the effective string tension the number of rephcas 
is expected to be finite. Indeed, the distance at which the string starts to break, L, will play 
the role of the infrared regulator in the formula similar to (jSHj), instead of fi. Therefore, as 
follows from Eq. (jHEj), ''^max ~ ^^{KqL) ~ 1, where we consider L ~ Kq^. In other words, 
confinement becomes less "binding" due to string breaking and the corresponding mass-gap 
equation supports fewer solutions — "bound states". A more detailed analysis of the QCD 
inspired interaction, including the proper string dynamics, from the point of view of replicas 
is in progress now and will be subject for future publications. 
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